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We report on the rapidity and centrality dependence of proton and anti-proton transverse mass 
distributions from 197 Au+ 197 Au collisions at */s~^ = 130 GeV as measured by the STAR experiment 
at RHIC. Our results are from the rapidity and transverse momentum range of \y\ < 0.5 and 0.35 
< pt <1.00 GeV/c. For both protons and anti-protons, transverse mass distributions become more 
convex from peripheral to central collisions demonstrating characteristics of collective expansion. 
The measured rapidity distributions and the mean transverse momenta versus rapidity are flat within 
\y\ < 0.5. Comparisons of our data with results from model calculations indicate that in order to 
obtain a consistent picture of the proton(anti-proton) yields and transverse mass distributions the 
possibility of pre-hadronic collective expansion may have to be taken into account. 



High energy nuclear collisions provide a unique op- 
portunity to study matter under extreme conditions for 
which one expects the formation of a system dominated 
by deconfined quarks and gluons lj. In the search for 



this deconfined state, baryons play an important role. 
Incoming beam baryons provide the energy for particle 
production and development of collective motion. It has 
systematically been observed that the net-baryon num.- 



3 



ber determines the chemical properties |2(. In addition, 
baryon transport and baryon production during the colli- 
sion are particularly interesting because of their dynam- 
ical nature 0, 0, IE IE 13 • However, these are difficult 
processes due to their non-perturbative features |9j, [lOj . 
At the RHIC energy ^/s^ =130 GeV, anti-proton to pro- 
ton ratios and yields at mid-rapidity have been reported 
by several experiments [llj, Il2l ll.'l Il4| . In the region of 
Pt ~ 2-3 GeV/c, the yield of protons approaches that 
of pions 01 m central collisions. The exact origin of 
this behavior is not clear and systematic measurements 
of baryon distributions are important. 

In this Letter, we present a systematic measurement 
of proton and anti-proton production in Au+Au colli- 
sions at y/SNN = 130 GeV in the rapidity range — 0.5 < 
y < 0.5 and for transverse momenta 0.35 < p t < 1.00 
GeV/c. In particular, we report the first RHIC mea- 
surements of the rapidity dependence of the proton and 
anti-proton yields, essential for exploring the existence of 
a boost-invariant region in the system. We also study the 
centrality dependence of the yields and mean transverse 
momenta for protons and anti-protons. These results al- 
low for a detailed comparison to model predictions of 
proton and anti-proton production at RHIC. 

Two independent 197 Au beams with an energy of 
65 GeV per nucleon were provided by the Relativistic 
Heavy Ion Collider (RHIC) at the Brookhaven National 
Laboratory. These beams collided around the geometric 
center of the Solenoid Tracker at RHIC (STAR). Charged 
particles stemming from these collisions were measured 
in a large volume Time Projection Chamber (TPC) [lfij . 
A large solenoidal magnet of 0.25 T field strength pro- 
vided momentum dispersion in the direction transverse 
to the beam line. 

For this analysis, we used 320k events with a minimum 
bias trigger and 154k events with a trigger selecting the 
10% most central events Events with a primary ver- 
tex within ±30 cm of the geometric center of the TPC 
along the beam axis were accepted. Tracks were required 
to have at least 23 out of 45 maximum possible space 
points in the TPC and to extrapolate back to the pri- 
mary vertex within 2 cm (distance of closest approach, 
dca) . To define the collision centrality, the measured raw 
multiplicity distribution of charged particles within the 
pseudorapidity range \rj\ < 0.75 was divided into eight 
bins. The highest centrality bin corresponds to 6% of the 
measured cross section for 197 Au+ 197 Au collisions (l6| . 
Protons and anti-protons were identified by correlating 
their energy loss dE/dx due to ionization in the TPC 
gas with the measured momentum. This method has al- 
ready been presented in [Tl| . 

The track reconstruction efficiency was determined 
by embedding simulated tracks into real events at the 
raw data level and subsequently applying the full re- 
construction algorithm to those events. The propaga- 
tion of single tracks was performed using the GEANT 
Monte Carlo code with a detailed model of the STAR 
geometry and a realistic simulation of the TPC response. 



The resulting track reconstruction efficiency is > 70% at 
Pt > 0.5 GeV/c for all centralities. By varying the track 
cuts, the overall systematic uncertainty in the track re- 
construction efficiency is estimated to be less than 10%. 
Further, the relative resolution in transverse momentum 
was derived to be w 4% at pt = 0.5 GeV/c. 

Secondary interactions of particles with the detector 
material generated background protons. Due to their 
different geometric origin, these background protons ap- 
pear as a rather flat tail in the cfca-distribution which 
extends into the peak region of primary protons at small 
dca. In order to correct for background protons, the pro- 
ton dca-distribution was fitted by the scaled anti-proton 
dca-distribution (which is background free) plus the re- 
sults on the proton background from Monte Carlo cal- 
culations. Raw yields were extracted for protons and 
anti- protons with dca < 2.0 cm, optimizing the signal to 
background ratio for protons. The raw yields were then 
corrected for track reconstruction efficiency, proton back- 
ground and in the case of anti-protons, for absorption 
in the detector material. The detector acceptance for 
protons(anti-protons) from the decay of lambdas(anti- 
lambdas) or other hyperons(anti-hyperons) is estimated 
to be larger than 95%. Corrections for feeddown from 
decays of hyperons(anti-hyperons) were not applied. 

The mid-rapidity (\y\ < 0.5) proton and anti-proton 
transverse mass distributions for all 8 centrality bins are 
shown in Fig. ^ Here, the transverse mass mt is given 

by m t — \Jp1 + TO p j with m p the rest mass of the pro- 
ton. The uncorrelated bin-to-bin systematic errors are 
estimated to be less than 7%. It is evident that both 
proton (left panel) and anti-proton (right panel) distri- 
butions become more convex from peripheral to central 
collisions indicating an increase in transverse radial flow. 
In order to extract pt-integrated yields, dN/dy and mean 
transverse momenta (pt), hydrodynamically motivated 
fits P"t| were applied, assuming a thermal source plus 
transverse radial flow. The fit parameters are the temper- 
ature Tf D at kinetic freeze-out and the transverse radial 
flow velocity j3 s at the system surface. A velocity profile 
Pt{r) — [3 s (r/R) - 5 was used, where R is the radius of the 
source. These fits simultaneously describe experimental 
spectra of charged pions 0]i kaons [x^, protons and 
anti-protons. The fit-results are shown as dashed lines 
in Fig. ^ The description of the experimental data is 
remarkably good. When strong collective flow develops, 
the transverse mass distributions for heavy mass particles 
will not have the simple exponential shape at low trans- 
verse mass. Therefore, the hydrodynamically motivated 
two parameter fits become necessary 20]. The increase 
of (pt) with centrality is indeed reflected in the values 
of the collective velocity parameter (Pt), which increase 
from about (0.42±0.10)c to (0.56±0.05)c from the most 
peripheral to the most central collisions, respectively. 

Note that in |lll |. the anti-proton transverse momen- 
tum distributions were fitted with a Gaussian function 
in p t . The difference between the model fit results and 
Gaussian fits in p t are < 6% and < 10% for (p t ) and inte- 
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grated yields dN/dy, respectively. Using other functions, 
i.e. exponential in m t and a Boltzmann function, the 
systematic uncertainty in dN/dy due to extrapolation is 
estimated to be less than 20%. Similarly, the systematic 
uncertainty in (p t ) is less than 6%. The total system- 
atic uncertainty in dN/dy is less than 22%, adding the 
contributions due to extrapolation (20%) and the track 
reconstruction efficiency (10%) in quadrature. The pro- 
ton and anti-proton rapidity distributions are shown in 
Fig. [21 (a) and (b) for different collision centralities. In 
the p t -range not covered by this experiment, the yield 
was extracted from the thermal model fit. The results 
are shown in Table I, which indicates that about 50% 
of the integrated yield was measured within the STAR 
TPC acceptance. The bin-to-bin systematic errors, due 
to background subtraction and PID contamination, are 
included in the plot. Since the shapes of the transverse 
mass distributions of protons and anti-protons do not dif- 
fer within statistical errors, the extracted values of (pt) 
shown in Fig. ^c) are the average of the two. Within 
| y | < 0.5, both values of (j>t) and dN/dy are found to be 
uniform as a function of rapidity indicating that at RHIC 
- for the first time in heavy ion collisions - a boost invari- 
ant region of at least one unit of rapidity for all centrality 
bins has developed. An analysis of charged particle ra- 
tios [2l| has demonstrated that at RHIC energies a boost 
invariant region does not exist at \y\ > 1.5. It will be of 
interest to study the rapidity distributions of different 
mass hadrons at RHIC. 

The top panels of Figure |31 show the (pt) within 
\y\ < 0.5 for protons (left) and anti-protons (right). The 
corresponding yields, dN/dy are shown in the bottom 
panels. The open symbols represent fiducial yields and 
filled ones show the integrated yields. The shaded bands 
indicate the systematic uncertainties in extracting (p t ) 
and dN/dy. Both values of (pt) and dN/dy are in good 
agreement with results from PHENIX Experimen- 
tal results on the lambda(anti-lambda) yields |23, show 
that the contribution of fccddown from hyperon decays 
to the proton(anti-proton) yields is «40%. The increase 
of (p t ) vs. centrality in the figure indicates the devel- 
opment of stronger collective expansion in more central 
collisions. Results from calculations with RQMD [23l |. 
RQMD with re-scattering switched off (w/o) and HI- 
JING pil I25I ] are represented by solid, dashed, and 
dashed-dotted lines, respectively. In the RQMD model 
[23L l2t| hadronic re-scattering has been implemented. 
This leads to the agreement with measurements in the 
mean transverse momentum. On the other hand, with- 
out the re-scattering, the HIJING model under-predicts 
the proton and anti-proton (pt), especially for central 
collisions. Overall, the model calculations fail to pre- 
dict the experimental yields consistently throughout the 
whole centrality range. Discrepancies between measured 
p/p ratios and predictions from RQMD and HIJING have 
been reported by other experiments [l3L Il4| . 

The bottom panels of Fig. 3 show that the ob- 
served mid-rapidity (\y\ < 0.5) proton and anti-proton 



yields, dN/dy are proportional to the number of charged 
hadrons. RQMD fails to predict the centrality depen- 
dence of the anti-proton yield due to the strong anni- 
hilation in hadronic re-scattering, especially in central 
collisions. Because of the annihilation, RQMD predicts 
a change in the p/p ratio of almost a factor of two from 
peripheral to central collisions, which is not consistent 
with observations [TT| . 

The results from RQMD reflect that within that model 
there is strong annihilation among baryons, and that 
large values of (pt) are built up from late hadronic rescat- 
terings. Based on RQMD, the annihilation of anti- 
protons created initially is expected to increase from 20% 
for perpiheral collisions, to 50% for the most central col- 
lisions. This is not consistent with the trend in Fig. 3, 
which indicates the measured proton and anti-proton 
yields increase approximately linearly with the number 
of negatively charged hadrons. This raises an important 
question. If, on the one hand the increase in annihila- 
tion with centrality predicted by RQMD is correct, then 
the centrality dependence of the initial baryon produc- 
tion must be much stronger than the linear dependence 
observed in Fig. 3, and the rough agreement between 
RQMD and the data for anti-protons is fortuitous. If, on 
the other hand, the agreement between RQMD and the 
linear dependence observed in Fig. 3 for anti-protons is 
correct, a possible explanation is that the anti-proton loss 
due to annihilation is smaller in central collisions than in 
peripheral collisions. This suggests the anti-protons may 
decouple from the surrounding matter early, and that the 
large experimental values of (pt) which are observed must 
arise from collective flow in the early stage [27], IH, |2t| . 
In order to distinguish this possibility from other possi- 
ble scenarios [3(J and study possible early-stage partonic 
collectivity at RHIC, systematic measurements of multi- 
strange baryons, charmed mesons, and particle correla- 
tions are necessary. 

In summary, we have reported on the centrality 
dependence of proton and anti-proton transverse mass 
and rapidity distributions from 197 Au+ 197 Au colli- 
sions at y/s NN = 130 GeV as measured by the STAR 
experiment at RHIC. The results reported here are 
from the rapidity and transverse momentum range of 
\y\ < 0.5 and 0.35< p t <1.00 GeV/c. For both protons 
and anti-protons, the transverse mass distributions 
become more convex from peripheral to central collisions 
indicating the enhancement of collective expansion in 
more central collisions. The rapidity distributions and 
(pt) versus rapidity are found to be flat within \y\ < 0.5 
suggesting a boost invariant region around mid-rapidity. 
The comparison of our data to results from microscopic 
transport models suggests that the observed collective 
expansion might have been dominantly developed at the 
early stage of the collision. 

We thank Drs. W. Busza, M. Gyulassy and V. Topor- 
Pop for exciting discussions. We wish to thank the RHIC 
Operations Group at Brookhaven National Laboratory 
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TABLE I: Mid-rapidity (\y\ < 0.5) proton and anti-proton 
results on (pt) and rapidity densities. The fiducial yield is 
measured within 0.35 < pt < 1.00 GeV/c. The errors are 
statistical. See text for discussions of systematic errors. 
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FIG. 1: Mid-rapidity (\y\ < 0.5) proton (left column) and 
anti-proton (right column) transverse mass distributions for 
most peripheral (bottom) to most central (top) collisions. The 
definitions of the centrality bins are listed in Table I. Rela- 
tively large systematic errors for protons in the low m t region 
are due to the background subtraction. Results from model 
fits are shown as dashed lines. 
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FIG. 2: The rapidity distributions of protons (a) and anti- 
protons (b) and the average transverse momentum {pt} (c), 
for most peripheral (bottom) to most central (top) collisions. 
The bin-to-bin systematic errors due to PID contamination, 
were included in the plot. Overall systematic errors due to 
extrapolation into the p t -range not covered by the experiment 
and the uncertainty in the track reconstruction efficiency are 
not shown in the figure. 
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FIG. 3: Mid-rapidity (pt) and dN/dy of protons and anti- 
protons as functions of the number of negatively charged 
hadrons. Open symbols are fiducial yields and filled ones are 
integrated yields. Systematic errors in the integrated yields 
are shown as shaded areas. Results from RQMD, RQMD with 
re-scattering switched off (w/o) and HIJING are shown as 
solid-lines, dashed-lines and dashed-dotted lines respectively. 
The experimental data and the results from RQMD and HI- 
JING include feeddown from hyperon decay. 



